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Description 

[0001] The invention relates to apparatus for measur- 
ing magnetic field infomiatlon, and maybe used in par- 
ticular for detennining components of magnetic field 
gradient. The system may be used to detect magnetic 
field gradients of the order of 1 00 fT m""" from a moving 
platfomri in the bacl<ground of the earth's magnetic field 
(approximately 70 \iT). 

[0002] SQUID (superconducting quantum interfer- 
ence devices) magnetometers are extremely sensitive 
devices which can measure vector components of mag- 
netic field as small as 1 0'^ times the ambient field of the 
earth. Measurement of fields near to the threshold of 
sensitivity in the presence of the earth's field present 
many difficulties. Geomagnetic noise and man-made 
noise, for example, are always dominant. Furthermore, 
unless the SQL) IDs are rigidly mounted, even minute 
motions In the earth's field will be reflected as gross 
changes in output and therefore It is often more useful 
to sense the gradient in a magnetic field, rather than the 
field itself. 

[0003] The gradient of a magnetic field may be meas- 
ured by using an intrinsic SQUID gradlometer An intrin- 
sic SQUID gradlometer has sensing colls made up of 
two loops connected with opposite polarity. In such a 
configuration the sensing coils must be highly balanced 
and aligned and small fractional changes in the effective 
size or orientation of each coil produce output signals 
from uniform fields which are indistinguishable from real 
gradients. Balancing is usually achieved by adjustment 
In special calibration rigs under laboratory conditions 
and is an expensive and time consuming process. Fur- 
thermore, as a large background field is incident on the 
sensor, currents are induced In the structure which give 
rise to 1/f noise and can give rise hysteresis problems. 
[0004] The gradient signal is larger for coils separated 
by a longer baseline, the useable baseline being limited 
by the tolerable inductance In the connections to the 
coils. This restriction on baseline can be removed by re- 
placing the gradlometer sensor by a pair of magnetom- 
eters, where outputs are subtracted to form a configured 
gradlometer. Such a configuration requires great stabil- 
ity and linearity in both time and frequency domains. 
Furthermore, as well as the difficulties in balancing, 
each sensor requires a very large dynamic range (better 
than 1 part in lO^) if it Is to be operable on a moving 
platform, as Is often required. In addition, there is still 
the problem that the field is incident on the structure. 
[0005] In US patent 5 1 22 744, an aligned three sen- 
sor configured gradlometer (Three SQUID Gradlometer, 
TSG) Is described in which a central sensor is used to 
feedback that component of the earth's field to colls sur- 
rounding each of the three sensors. The outputs of the 
outer pair of sensors are subtracted and this difference 
gives a measure of the required gradient if the feedback 
field is uniform. In this configuration, the dynamic range 
is considerably Improved and one component of the 



earth's field is not directly incident upon the sensors. 
[0006] The same technique has been extended to 
gradiometers based on fluxgate magnetometers rather 
than SQUID magnetometers [ R,H. Koch et at: "Room 

5 temperature three sensor magnetic field gradiometer" 
Review of Scientific instruments, Jan. 1996, AIR USA, 
vol. 67, No. 1, pages 230-235]. 
[0007] However, the problems of stability, non-linear- 
ity, unifomnity of feedback and the need to calibrate and 

10 fix the balance are not overcome. It is the difficulty of 
accurately subtracting the sensor outputs which give 
rise to many of the problems. 

[0008] Also relevant to the present invention is the 
processing technique used to compensate for motion 
15 noise in extremely low frequency submarine receiving 
aniennaslR.J.Dingerand J.R.Davis, Proc. IEEE, vol.64, 
No. 10, Oct. 1976\. 

[0009] An additional technique for improving the per- 
fonnance of SQUID sensors, based on adaptive positive 

20 feedback, is known from US patent 5 488 295. 

[0010] For the purposes of this specification, thetenn 
magnetometer shall be taken to refer to a device for 
measuring the magnetic field component in a particular 
direction and the term gradiometer shall be taken to re- 

25 fer to a device for measuring magnetic field gradient 
components. A total field magnetometer shall be taken 
to refer a device for measuring the total magnetic field 
i.e. the square root of the sum of the squares of the mag- 
netic field components in three orthogonal directions. 

30 [0011] According to the present invention, a system 
for measuring one or more magnetic field gradient com- 
ponent of a magnetk: field comprises; 

(i) at least two magnetic sensors for sensing a mag- 
55 netic field, wherein each sensor generates a sensor 

output, said sensor outputs having an associated 
total energy, E, and wherein at least two of the sen- 
sors are arranged to sense the magnetic field In 
substantially the same direction and 

40 

(ii) means for performing adaptive signal process- 
ing of the sensor outputs such that the system is 
adaptively balanced, whereby said means generate 
one or more magnetic field gradient components, 

45 

characterised in that 

(iil) the system incorporates global feedback means 
for providing a substantially uniform magnetic field 
50 at the two or magnetic sensors, and 

(iv) the adaptive signal processing means comprise 
means for minimising the total energy, E, of the sen- 
sor outputs subject to a constraint, whereby the 
55 constraint determines which of one or more mag- 
netic field gradient components is generated. 

[0012] The invention provides the advantage that it 
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eliminates the need to calibrate and fix the balance un- 
der controlled conditions and to maintain that calibration 
for long periods. Furthemriore, the requirements on the 
mechanical rigidity and stability required for low noise 
operation are considerably relaxed. The advantages of 
the known configured systems are maintained. 
[0013] In a preferred embodiment, the means for min- 
imising the total energy, E, of the sensor outputs may 
also generate a total magnetic field measurement, 
whereby the constraint determines which of one or more 
magnetic field gradient components or a total magnetic 
field measurement is generated. 
[0014] In a further preferred embodiment, the gradi- 
ometer may comprise a computer on which an adaptive 
signal processing algorithm (ASPA) is loaded. 
[0015] The gradiometer global feedback means may 
also be arranged to provide a substantially unifonn mag- 
netic field gradient at the two or more magnetic sensors. 
The global feedback means may comprise at least one 
global feedback coil set. For example, each set may 
comprise two or more Helmholz coils. 
[0016] The gradiometer may also comprise; 

means for generating at least one difference signal 
between two sensor outputs, wherein said sensor 
outputs each correspond to a magnetic field In sub- 
stantially the same direction, and 

analogue to digital conversion means for converting 
the one or more difference signal and the two or 
more sensor outputs into equivalent digital data. 

[0017] At least one of the magnetic sensors may be 
one of a fluxgate, a Hall probe, a magneto-resistive sen- 
sor or a superconducting quantum interference device 
(SQUID) magnetometer. Alternatively, at least one of 
the magnetic sensors may be a gradiometer. 
[0018] If a SQUID magnetometer is included in the 
gradiometer, the gradiometer also includes cooling 
means for reducing the temperature of the SQUID mag- 
netometer. Each SQUID magnetometer may have as- 
sociated local feedback means for maintaining a sub- 
stantially constant state of magnetic flux in the respec- 
tive SQUID magnetometer 

[0019] In one embodiment of the invention, the gradi- 
ometer comprises; 

at least four magnetic sensors for sensing a mag- 
netic field, wherein three of the sensors are ar- 
ranged such that they sense the magnetic field In 
three substantially orthogonal directions and 
wherein at least two of said sensors are arranged 
such that they sense the magnetic field in substan- 
tially the same direction. 

[0020] In a further preferred embodiment, the gradi- 
ometer may comprise at least eight magnetic sensors. 
For example, the eight or more magnetic sensors may 



be arranged at the vertices of a tetrahedron structure. 
In this configuration, the three global feedback coil sets 
may be oriented in three substantially orthogonal direc- 
tions, for generating a substantially unifonn magnetic 
5 field at each of the eight or more magnetic sensors. 
[0021] In another embodiment, the gradiometer may 
comprise at least three magnetic sensors, wherein at 
least three of the sensors are oriented in substantially 
the same direction and whereby the means for perform- 
10 ing adaptive signal processing may generate a magnet- 
ic field gradient component of at least second order. 
[0022] in another embodiment of the invention, the 
system may be arranged to provide a biomagnetic sens- 
ing system. The at least one global feedback coil set 
15 may be an-anged to surround a subject, for example a 
human subject, generating a magnetic field to be meas- 
ured wherein the subject may be in close proximity to 
the magnetic sensors. This system provides an advan- 
tage over conventional biomagnetic sensing systems in 
20 that the large magnetically shielded room required in 
conventional systems is no longer needed. 
[0023] According to another aspect of the invention, 
a method for measuring at least one magnetic field gra- 
dient component using a gradiometer comprises the 
25 steps of; 

(i) sensing a magnetic field component at two or 
more positions using two or more magnetic sen- 
sors, wherein at least two of the magnetic field com- 

30 ponents are sensed in substantially the same direc- 
tion, 

(ii) generating two or more output signals, having 
an associated total energy, E, corresponding to said 

35 magnetic field components, 

(iii) providing global feedback in the fonm of a sub- 
stantially unifonnfi magnetic field at the two or more 
magnetic sensors, 

40 

(iv) performing adaptive signal processing of the 
output signals and minimising the total energy, E, 
of the output signals subject to a constraint such 
that the gradiometer is adaptively balanced, 

45 

(v) constraining the minimisation of the total energy, 
E, such that a magnetic field gradient component is 
generated, and 

so (vi) generating at least one. magnetic field gradient 
component measurement. 

[0024] In one embodiment of this aspect of the inven- 
tion, the method may comprise the steps of; 

55 

(i) minimising the total energy, E, of the output sig- 
nals subject to a constraint such that the gradiom- 
eter is adaptively balanced, and 
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(ii) constraining the minimisation of the total energy, 
E, such that a total magnetic field measurement is 
generated. 

[0025] In an alternative embodiment of this aspect of 
the invention, the method may comprise the step of pro- 
viding a substantially uniform magnetic field gradient at 
the two or more magnetic sensors. The method may be 
used for measuring magnetic field components in a hu- 
man subject in close proximity with the two or more mag- 
netic sensors. 

[0026] In another embodiment of this aspect of the In- 
vention, the method may comprise the steps of; 

(i) sensing the magnetic field component at three or 
more positions, wherein at least three of said mag- 
netic field components are sensed in substantially 
the same direction, 

(ii) minimising the total energy, E, of the output sig- 
nals subject to a constraint such that the gradiom- 
eter is adaptively balanced, and 

(iii) constraining the minimisation of the total energy, 
E, such that a magnetic field gradient component of 
at least second order is generated. 

[0027] The system utilises an adaptive signal 
processing technique to achieve high levels of balance. 
The application of this technique to magnetic detection 
with gradiometer systems is unl<nown. The employment 
of adaptive signal processing is essential to allow the 
system to utilise a scheme of global (or overall) magnet- 
ic field feedback. For a system from a moving platform, 
this provides sufficient dynamic range for the system to 
operate In the earth's field whilst maintaining high levels 
of balance. Furthermore, it also provides the advantage 
that it prohibits the ambient external field impinging di- 
rectly on the sensors. Furthermore, it relaxes the me- 
chanical rigidity constraints for the system. 
[0028] The invention will now be described, by exam- 
ple only, with reference to the following figures in which; 

Figure 1 shows a schematic diagram of a conven- 
tional configured gradiometer device for measuring 
magnetic field gradients, based on two SQUID mag- 
netometers, 

Figure 2 shows a schematic diagram of a three 
SQUID gradiometer device known in the prior art, 

Figure 3 shows a diagram of a SQUID magnetom- 
eter device, 

Figure 4 shows a flow diagram to illustrate an adap- 
tive signal processing algorithm (ASPA) which may 
be included in the invention, 



6 

Figure 5 shows a schematb diagram of a single axis 
global feedback arrangement of the present inven- 
tion, 

5 Figure 6 shows a schematic of a four sensor em- 
bodiment of the present invention to illustrate a 
means for obtaining good balance to measure one 
gradient, 

10 Figure 7 shows the results obtained from an oper- 
ational system arranged as in Figure 6. 

Figure 8 shows a schematic diagram of a global 
feedback gradiometer for obtaining five independ- 
is ent magnetic field gradient measurements and 

Figure 9 shows an arrangement of a biomagnetic 
sensing system of the invention. 

20 [0029] Refemng to Figure 1, a conventional config- 
ured gradiometer for measuring the magnetic field gra- 
dient may comprise two magnetometers 1 a, 1 b, typically 
SQUID magnetometers, separated by a distance d 
(baseline distance), whereby each magnetometer 1 a,1 b 

25 measures the magnetic field at its particular location. 
The gradiometer also comprises amplifiers 2a,2b and 
integrators 3a,3b, feedback coils 4a,4b and resistors 5a, 
5b, having the same resistance, for providing a correc- 
tion current producing a field equal and opposite to that 

30 of the external magnetic field incident on the SQUIDs. 
Voltages Va, Vb are output from the device correspond- 
ing to the magnetic field at each of the SQUID magne- 
tometers 1 a,1 b respectively. The difference between the 
two voltages, Va-Vb, provides an approximation of the 

35 magnetic field gradient. 

[0030] The two magnetometer gradiometer shown in 
Figure 1 has a large common mode signal of the two 
SQUIDs, arising from the earth's magnetic field and re- 
quires an almost impossibly large degree of common 

40 mode rejection of an amplifier taking the difference In 
the outputs of the two SQUID magnetometers 1a, 1b. 
Furthermore, each magnetometer and associated am- 
plifier requires a very large dynamic range (better than 
1 part in 10^) and an exceptionally high linearity if it is 

45 to be successfully operated from a moving platfomi. 
[0031] Refemng to Figure 2, a conventional three 
magnetometer configured gradiometer (TSG) (US pat- 
ent 5 122 744) comprises three SQUID magnetometers 
6,7,8 and five feedback coils 9,10,11,12,13. The 

50 SQUIDs 6,7 are sensor SQUIDs used to measure the 
magnetic field at their respective locations and are sep- 
arated by a distance /. The SQUID 8 acts as a reference 
magnetometer and is used to cancel out the effect of the 
background magnetic field from the earth's magnetic 

55 field In SQUIDs 6,7 via the feedback coils 9,10 provided 
with a signal from the reference SQUID 8. Amplifiers 
1 4,1 5,1 6, and integrators 1 7,1 8, 1 9 are also shown in the 
figure. 
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[0032] The reference SQUID magnetometer 8 is op- 
erated in a similar way to the sensor SQU ID magnetom- 
eters 6,7 except the feedback signal for SQUID 8 is ap- 
plied to all three SQUIDs 6.7,8 via the feedback coils 
9,10,11. Amplifier 6 and integrator 1 9 are used to gen- 5 
erate the feedback signals to coils 9,10,11 whteh sub- 
tract the background magnetic field from SQUIDs 6,7.8 
respectively. 

[0033] The field gradient is measured by subtracting 
the outputs from the two sensor SQUIDs 6,7. One ad- 
vantage of the three SQUID system is that these outputs 
are much smaller than those measured using a two 
SQUID gradiometer (Figure 1). However, in this config- 
ured system the problems of stability, non-linearrty and 
uniformity of feedback, and the need to calibrate and fix 
the balance, are not overcome. 

[0034] By way of background to the invention, the op- 
eration of a SQUID magnetometer is described with ref- 
erence to Figure 3. The SQUID 20 comprises a super- 
conducting loop 21 having at least one weak link 22 (e. 
g. Josephson junction) which can exhibit a periodic cur- 
rent-phase relation. The SQUID is magnetically coupled 
to a superconducting flux transformer via an input coil 
23. The flux transformer pickup coil 24 senses the field 
to be measured. When there Is a change in the magnetic 
field to be detected, a circulating current will be induced 
in the input coil 23. The circulating current produces a 
magnetic field which couples to the SQUID loop 21 and 
is detected. A SQUID (as shown in Figures 1 ,2,3 and 4) 
would usually have an associated flux transformer 
[0035] The gradiometer of the present invention can 
be configured to measure magnetic field gradients, and 
magnetic fields, originating from distant sources in the 
presence of the large and substantially uniform field of 
the earth, in particular if the gradiometer is operated 
from a moving platfomn. Alternatively, the gradiometer 
may be configured to measure magnetic field gradients 
or magnetk: fields originating from sources which are in 
close proximity to the sensor, where the sensor is not 
moving but is operating in a magnetically noisy environ- 
ment. 

[0036] Figure 4 shows one embodiment of the inven- 
tion; a single-axis global feedback gradiometer. The 
system comprises two magnetometers 25a.25b and 
may be used to measure either magnetic field gradients 
or, in some configurations, the total magnetic field. The 
magnetometers 25a.25b are represented as SQUID 
magnetometers, although other magnetometer devices 
may also be used, for example, flux gate, Hall probe 
sensors or magneto-resistive devices. 
[0037] The two magnetometers 25a, 26b are situated 
one above the other (axial configuration) or side by side 
(off-axis configuration) and sense the magnetic field in 
a particular direction. In the case of SQUID magnetom- 
eters, each SQUID 25a,25b has its own local feedback 
coil 26a,26b which surrounds the associated magne- 
tometer 25a,25b. These local feedback coils 26a,26b 
provide the flux locking for the respective SQUID mag- 



netometers 25a,25b, as known in the prior art (see Fig- 
ures 1 and 2). 

[0038] Each channel has an associated amprrfier27a, 
27b and integrator 28a,28b and therefore each magne- 
tometer 25a.25b has readout electronics which operate 
in a flux locked loop mode, using integral feedback via 
28a,28b, such that the electronics' feedback flux main- 
tains a constant state of flux, whether finite or zero, in 
each magnetometer 25a,25b. 

[0039] As the flux in the magnetometers 25a ,25b is 
constant, the feedback current is proportional to the ex- 
ternal field applied to each magnetometer 25a,25b and 
the outputs 29a,29b from the magnetometers' readout 
electronics therefore give a measure of the feedback 
current in the relevant magnetometer channel. The 
magnetometer channel outputs 29a, 29b, from integra- 
tors 28a,28b respectively, are then input to adaptive sig- 
nal processing means (to be described later). 
[0040] The system also comprises global (common 
mode) feedback coils 30a,30b whteh enclose the mag- 
netometers 25a,25b. For the purpose of this specifica- 
tion, the term global feedback shall be taken to mean 
the use of one or more remote coil sets to apply mag- 
netic feedback to all of the magnetic sensors simultane- 
ously. The output signals 29a,29b from the magnetom- 
eter channels are passed to a component 31 for elec- 
tronically adding the magnetic field outputs from the two 
magnetometers 25a,25b and then to a global feedback 
integrator 32. The output 33 from the global feedback 
integrator32 is passed through the global feedback coils 
30a,30b and may be passed to the signal processing 
means. 

[0041] The feedback field realistically emulates the 
real external field and is applied approximately equally 
to both magnetometers 25a, 25b. Enclosing the magne- 
tometers 25a,25b within the global coils 25a,25b there- 
fore enables the earth's magnetic field to be suppressed 
at the magnetometers 25a,25b. The output from the glo- 
bal feedback integrator 32 gives a measure of the cur- 
rent in the feedback coils 30a,30b. The output data, for 
input to the data processing means, therefore compris- 
es the individual magnetometer outputs 29a,29b and 
the current in the global feedback coils 30a,30b (output 
33). 

[0042] The use of global feedback eases the match- 
ing requirements of the separate magnetometer chan- 
nels and eliminates problems associated with field in- 
duced noise, hysteresis and field modulation problems 
are also reduced. It also eases the required orientational 
alignment accuracy of the magnetometers 25a,25b. The 
use of the local feedback coils 26a, 26b associated with 
each SQUID magnetometer 25a,25b enables the low 
noise and fine resolution required for magnetic field gra- 
dient measurement. Furthermore, it provides enough 
dynamic range in the global outer loop 30a,30b to han- 
dle the earth's magnetic field. 

[0043] When operated from a moving platfonm to 
measure distant magnetic sources, the magnetometers 
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25a,25b to be able to cope with being turned over in the 
earth's magnetic field (-70 \lT), which represents a 
change of -140 ^iT Typically, SQUID magnetometers 
are sensitive to 10 ff and therefore this is a dynamic 
range of 1 part in 1.4 x 10^°, or 203 dB, which is not 
achievable using semiconductor based readout elec- 
tronics. However, by cancelling out the unifonn field and 
measuring magnetic field gradients the dynamic range 
problem may be overcome. 

[0044] When operated from a fixed platform to meas- 
ure near magnetic sources, the cancellation of unifomi 
field is also preferable. In such measurements, distant 
magnetic gradient sources may also interfere with the 
measurements. This may be overcome by measuring 
and cancelling the extemal gradients using global feed- 
back. 

[0045] The following description relates to the gradl- 
ometer when configured to measure distant sources 
from a moving platfonn. Consider the two magnetome- 
ters 25a,25b sensing magnetic fields and Bq. The 
field gradient is (B;^ - Bg) and the average field half way 
between them is (B^ + 83) / 2 which is proportional to 
the sum of the fields. The magnetic fields B^, Bg are 
summed (31), and this sum is fed back, via Integral feed- 
back to the global coils 30a,30b (the factor 1/2 is adjust- 
ed for automatically by the integrator 32). The magnetic 
field half way between the magnetometers 25a,25b is 
therefore always maintained constant (or zero if initially 
zero), as required. 

[0046] The two magnetometers 25a,25b now have 
the uniform component of the earth's magnetic field, Bg, 
removed, to give (B^ - Bg) and (Bq - Bg). If the magnetic 
field gradient is required, (B^ - Bg), in conventional 
gradiometer systems this is derived by subtracting the 
magnetometer outputs. However, it is almost certain 
that the magnetometers will have slightly different scal- 
ing of their outputs (mismatch) and be very slightly mis- 
aligned and misoriented. as conventional electronics 
may only be linear to about 1 in 1 0^ and may be matched 
with similar accuracy. Misorientation of the magnetom- 
eters by only 0.1 degree for example leads to an error 
of 1 .5 X 1 0-6, In order to measure magnetic fields of the 
order of 1 GO fT m"'' in the background of the earth's mag- 
netic field (- 70 \iJ), an accuracy of the order of 1 part 
in 10^0 is required and this is clearly difficult to achieve 
with such misorientation effects. 
[0047] Rather than attempting to subtract the sensor 
outputs, as In the prior art, the present system utilises 
an energy minimisation algorithm which Is subject to a 
suitable constraint designed to recover the gradient In- 
formation. An energy minimisation algorithm may also 
be employed to obtain a total field measurement from 
three orthogonal magnetometers in the same way. 
[0048] This method of obtaining the gradient infomia- 
tion is adaptive In the sense that It constantly tries to 
improve the balance based on the history of the sensor 
outputs over a chosen timescale. The use of this energy 
minimisation algorithm allows the use of global (com- 



10 

mon mode) feedback to overcome many of the prob- 
lems in existing systems, such as the problems of sta- 
bility, unlfomilty of feedback and the need to calibrate 
and fix the balance. 

5 [0049] An adaptive signal processing algorithm (AS- 
PA) is used to obtain the magnetic field gradient meas- 
urement from the magnetometer outputs and is an es- 
sential aspect of the invention in order to achieve the 
high level of balance required. Furthennore, the AS PA 

10 is advantageous In that It offers a balancing technique 
which may be perfonned continuously in real time. The 
principle of the ASPA Is based on a well known mathe- 
matical techniques which is used extensively in other 
types of system, for example adaptive beam fonnlng in 

^5 phased array antennas. Further details of this technique 
may be found in Radar Array Processing, by S. Haykin, 
J. Litva and T J. Shepherd (pub. Spinger-Veriag 
(1993)), 

[0050] In the case of gradient measurement the ASPA 
functions by fomiing the outputs of (for example) two 
magnetometers, aligned and separated In the z direc- 
tion, B21 and B22 Into a linear combination, co^B2-i + 
(02822- The energy of this combination, over a defined 
time window, Is minimised by adjusting the coefficients 
(Di and (D2 subject to one or more constraints. The pur- 
pose of the constraint is to avoid the trivial solution 
co^=CB2=0 and yield the desired gradient information. 
This procedure is referred to as the adaptive balancing 
of the gradiometer system. 

[0051] The fonn of the constraint Is co = \l, where 
S and (1 are chosen constraint matrices ("'' = transpose). 
For example, if the two magnetometers 25a,25b were 
perfectly matched and aligned, and the following con- 
straint is applied; 



(where (1 0) and ^ = 1) i.e. co^ = 1 , the algorithm 
would calculate 0)2 = -1 . The output, which Is the linear 
combination <o^B^^ + 0)2622, would then be B21 - B22 i. 
45 e. the required gradient. If the sensors were mismatched 
then will automatically differ from -1 by that mis- 
match, thus con-ecting for it without the need for calibra- 
tion. 

[0052] A flow diagram to Illustrate the operation of the 
50 ASPA is shown in Figure 4. The output at time t, y(t), is 
given by; 

y(t)=x^(t)^ 

55 

\Nherex(t) Is the Input data set at time f (i.e. the magne- 
tometer outputs) and © is given by M(t)'^ S(S^M(ty^ S)-^ 
\i. M(t) is the covariance matrix at time t and is given by; 
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[0053] For each new input data set, x(t), the algorithm 
connputes the covariance matrix, M(t), which depends 
on the preceeding covariance matrix, MIU 1), at time M . 
[0054] The system shown tn Rgure 5 Includes two 
magnetometers 25a,25b for sensing the magnetic field 
and Illustrates the principle of operation of a gradiometer 
system with global feedback which may be used to give 
infomiation on both the field gradient in one direction 
and the total magnetic field e.g. 8^ and dB/dx. If the two 
magnetometers 25a,25b are perfectly aligned, two mag- 
netometers are sufficient to give an accurate field gra- 
dient measurement. However, if there is any misalign- 
ment of the magnetometers (which inevitably occurs), 
the magnetometers 25a,25b may also be sensitive to 
magnetic fields in the orthogonal y and z directions. 
[0055] I n practice, therefore, a m ore useful system will 
Include magnetometers for monitoring the two direc- 
tions orthogonal to the gradient sensing magnetometers 
25a,25b. These fields are then accounted for in the AS- 
PA to give a more accurate measure of the field gradient 
dBj/dx. (In addition, the measurements of By and 
may be used with that of B^ to obtain the total field using 
a separate AS PA.) 

[0056] This is illustrated in Figure 6 which shows the 
magnetic field components B21, B22, Bxi, Byi which are 
measured by four sensors (not shown in the diagram). 
The four measured components, B21, B22, Bxi, Byi are 
Input to an ASPA 34 (via an Analogue to Digital Con- 
verter (ADC)35). The inputs to the ASPA 34 are as- 
signed the weights (o^ , oDg, (o^, and 004 respectively. One 
constraint could be 



(10 0 




which fixes the coefficient co^ to 1 and allows all the oth- 
ers to become freely adjustable by the energy minimi- 
sation. The result would be that would adjust to neariy 
-1 . as described eariier, and CO3 and CD4 would become 
the required small numbers which correct for misalign- 
ment. The output of the ASPA 34 Is then; 
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i.e. the required gradient. 

[0057] An alternative method which may be prefera- 
ble would be to perform a subtraction of B^i and 3^2 
and use that as one input to the ASPA 34, The input 
matrix would then be (821*622, Byv Bz^) and the 
same constraint would apply. 

[0058] The energy minimisation occurs over a time 
frame, set in the software, which is long compared with 
the signal to be measured. The system therefore be- 
comes self balancing, adjusting itself over that time 
frame. Any changes of balance comparable with or 
slower than that period are adapted out, therefore relax- 
Ing the stability constraints which apply in simple sub- 
traction methods. Any magnetic field gradients chang- 
ing more rapidly than the adaptation time will be appar- 
ent as measured gradient signals. The algorithm has an 
efficient update procedure such that each new set of 
points allows the window to slide along the time series. 
[0059] The system shown in Figure 6 requires four 
magnetometers to accurately measure one gradient, al- 
though for each additional gradient measurement only 
one further magnetometer is required. For example, the 
addition of one further magnetometer enables a meas- 
urement of dBy/dy, with non-idealities in this gradient 
corrected for by measurements of 8^ and B^ 
[0060] The potentially most useful configuration of the 
system is a 5-axis gradient gradiometer system. This 
system may be used to measure the minimum number 
of gradients required to model all the independent gra- 
dients of the magnetic field due to a dipolar source. As 
a minimum requirement, a 5-axrs gradient gradiometer 
system for measuring gradients dB/dx, dB/dy, dB/dy, 
dBj/dz and dBy/dz comprises eight magnetometers in 
total. In practice, however, it may be preferable to allow 
for some redundancy and to configure a 5-axls gradient 
system with more than the minimum requirement of 
eight magnetometers. 

[0061] A multiple-axis gradiometer system operates 
on the same principle as the single-axis gradiometer de- 
scribed previously. Typically, the global feedback may 
be applied through a coil set comprising three main 
feedback coils, such as Helmholz coils or higher order, 
one for each direction, x,y.z. The coils set may also com- 
prise additional minor trimming coils to improve the 
matching at the magnetic sensors and to overcome any 
nonunifomriity In the global colls. The outputs to the coils 
may be used as the orthogonal field estimators for a total 
tield evaluation by an ASPA. A similar procedure is used 
to determine the magnetic field gradient measurements. 
A separate but simultaneous energy minimisation is per- 
formed for each gradient component. The input data are 
either the time series from each magnetometer or linear 
combinations of the magnetometer outputs. The con- 
straints are then chosen to force the choice of different 
gradiometer components. 

[0062] Figure 7 shows the results obtained for an op- 
erational system, as illustrated schematically in Figure 
6, comprising two SQUIOs in the Z direction (Z^ and Z2 
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respectively) and one in each of the orthogonal direc- 
tions (X and Y). The upper four traces show outputs B^^, 
B22, Bxi, Byi {as in Figure 6) and the lower trace shows 
the gradient measurement B21-B22. The outputs are 
generated by the movement of a small magnetic dipole. 5 
The opposing nature of the and Z2 signals is an indi- 
cation that the global feedback scheme Is operating cor- 
rectly i.e. the field half way between the Zi and sen- 
sors is nulled. 

[0063] Figure 8 shows a schematic diagram to illus- 10 
trate how a g radiometer system may be configured to 
provide a total field output measurement, B^, and five 
gradient outputs. G^z (dB/d^, G^^f {dB/d)^, (dB/ 
d!z), Gyx (c/By/cfx) and G^ (dB/dz). In this example, the 
gradiometer comprises nine magnetometers (not 15 
shown for clarity), the outputs from which provide meas- 
urements of the components of magnetic field, B^^ 0^ 
Sjfj, Byp By2, Byj, B^p B^ B22 (the directions of the 
components of magnetic field are shown in the figure). 
The magnetic field measurements in each direction, x, 20 
y and z, are summed electronically and are passed to 
associated integrating circuitry 36a,36b,36c, the out- 
puts from which pass through transductance amplifiers 
37a, 37b, 37c and provide a current to three global feed- 
back coils 30a,30b,30c respectively. 25 
[0064] The sums of the magnetic field In each direc- 
tion x,y,z (i.e. the outputs from the integrators 36a,36b, 
36c) give a measure of the total magnetic field in the x, 
y and z directions, B^^ By, and B^ respectively. The glo- 
bal feedback coils 30a,30b,30c are oriented in substan- 30 
tially orthogonal x,y and z directions and, as described 
previously for the single-axis gradiometer system, have 
the effect of cancelling the ambient magnetic field. 
[0065] In order to derive the five gradient measure- 
ments, G^^, Gjof, G^^, G^ and G^, the magnetic field 35 
gradient measurements, B^^-B^, B^-B^ ^xi'^x2' ^y2' 
Sy^and By2-Byp are deduced electronically (not shown 
in the figure) from the magnetometer outputs and these 
gradient estimators are passed to the ADC 21. The 
equivalent digital outputs 42,43,44,45,46 from the ADC 40 
35 are passed to an ASPA 34b, loaded on a computer, 
along with the digital outputs 39,40,41 , for adaptive sig- 
nal processing in order to generate the accurate mag- 
netic field gradient information. 

[0066] In order to obtain a measurement of the total 45 
magnetic field, By, the integrated outputs. B^ By, and 
B^ are passed to an ADC 35 and the equivalent digital 
outputs from which 39,40,41 are passed to an adaptive 
signal processing algorithm (ASPA) 34a. 
[0067] For clarity, two adaptive signal processing al- so 
gorithms 34a,34b are shown in Figure 8. In practice, 
however, a separate energy minimisation algorithm 
(ASPA) is required for each independent gradient meas- 
urement. For example, to detennine G^^, the sums of 
the magnetic field components in each direction, x, y 55 
and z (digital signals 39,40,41 ) and the difference signal 
^zr^z2 (digital signal 42) are used for energy minimisa- 
tion. A separate energy minimisation routine is also re- 



quired for the total field measurement. In practice, how- 
ever, it is be convenient to load all the required algo- 
rithms onto a single computer 
[0068] The example shown in Figure 8 represents on- 
ly one way of linearly combining the outputs from the 
nine magnetometers so as to obtain five magnetic field 
gradient measurements. Other linear combinations of 
the magnetometer outputs may also be used as inputs 
to the adaptive signal processing algorithm to obtain 
such infomnation. 

[0069] The rigidity requirements on the complete 
structure mean that preferably it should be as light and 
stiff as possible. One suitable construction may be a tet- 
rahedron of thin walled carbon fibre tubes with groups 
of magnetometers placed at each vertex. Alternatively, 
a cruciform structure may be used. 
[0070] For a gradiometer system designed to be 
moveable, e.g. on a moving platfomn, the baseline may 
typically be between 40 mm and 100 mm, as smaller 
structures give better rigidity than larger. However, if the 
gradiometer is to remain static, a larger baseline may 
be used to give a less rigid gradiometer with better sen- 
sitivity. The same sensitivity may be achieved by using 
more sensitive magnetometers on a shorter baseline to 
increase rigidity, for example SQUID magnetometers 
fabricated on 25 mm chips at 40 mm spacing, rather 
than 10 mm chips on a 100 mm baseline. 
[0071] Preferably, the global feedback coils must be 
large enough to generate a substantially unifonn field at 
the magnetometers and the size of the global feedback 
coils is determined to some extent by the particular ap- 
plication for which the system is to be used. For exam- 
ple, if a large space is available, perhaps for a static 
structure where the baseline is greater, the global feed- 
back coils may be large, typically 1 m In diameter, cre- 
ating a uniform field at all magnetometers. Conversely, 
if the baseline is small and the available space is limited, 
the global feedback coils may be much smaller, for ex- 
ample 1 0-20 cm. Particulariy for global coils of smaller 
diameter, minor trimming colls may be used to Improve 
the field matching at the magnetometers and overcome 
any non unifomnity in the global coils. 
[0072] The importance of having global feedback in 
the system depends on the sensitivity of the particular 
magnetometers Included in the system. The more sen- 
sitive the magnetic sensor, the more important it is to 
use global feedback to provide a high dynamic range. 
For the system shown in Figure 4, comprising SQUID 
magnetometers of high sensitivity, it is preferable to in- 
clude global feedback in the system. However, other 
magnetic sensors may be included in the system, for 
example flux gates, Hall probes or magneto-resistive 
devices, which are typically less sensitive than SQUID 
magnetometers but may be sufficiently sensitive for cer- 
tain applications. 

[0073] If less sensitive magnetometers are included, 
and there is a lower requirement on the dynamic range, 
the system may be operated without global feedback. 
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In this case, the magnetometer outputs are input to the 
ASPA for adaptive signal balancing, as described pre- 
viously, where the required magnetic field gradient or 
total field infomnation may be output. 
[0074] If SQUID magnetometers are included in the 
system they must be operated within a cryogenic envi- 
ronment (e.g. a cryostat containing liquid nitrogen or he- 
lium) to maintain the SQUIDs at a suitably low temper- 
ature. The local feedback coiis may be Integrated on the 
SQUID magnetometer chips which typically have a di- 
mension of between 1 0 mm and 50 mm. The local feed- 
back coils provide the flux locking for the SQUID. The 
use of two sets of nested feedback loops achieves the 
requirement of low noise and fine resolution for the pur- 
poses of measuring field gradient, while providing suffi- 
cient dynamic range in the outer global feedback loop 
to handle the earth's magnetic field. 
[0075] In an altemative embodiment, the system may 
be used to measure second or higher order magnetic 
field gradient components (i.e. <fB^d)^. <fB^<i^e\c.), 
For example, the two magnetometers shown in Figure 
5 may be replaced with three magnetometers orientated 
in substantially the same direction and spaced equidis- 
tantly to measure B^.,, and Bj^g. The three outputs 
could be summed and integrated to provide the required 
global feedback. The required linear combination to 
measure d^B^/dx^ is Bxi-2Bx2*Bx3 which is easily ob- 
tained by choosing a suitable constraint in the ASPA. 
This configuration may be extended to a multiple axis 
system. In particular, a system for measuring second or 
higher order magnetic field gradients may be useful for 
applications requiring exceptionally high rejection of ex- 
ternal interference, such as medical scanning applica- 
tions. 

[0076] For some applications, the global feedback 
coils may be substantially larger than the magnetometer 
array. This configuration would be useful where the 
source of the gradient fields to be measured is in close 
proximity to the magnetometer array. An example of this 
is in a biomagnetic sensor system. For example, refer- 
ring to Figure 9, a human subject 50 in a room 55 may 
be examined for magnetism using an array of magne- 
tometers 25, such as an array of SQUID magnetome- 
ters. In the arrangement shown, global feedback coils 
30a,30b,30c surround the human subject 50 and the 
magnetometers 25. The global feedback coils 30a. 30b, 
30c are arranged in three substantially orthogonal 
planes and provide a substantially uniform magnetic 
field atthe magnetic sensors 25. In addition, global feed- 
back coil sets 60a.60b,60c may be configured to feed 
back the first order magnetic field gradients (the front 
most component of coil set 60c is not shown for clarity 
and part of the rear most component is obscured by the 
human subject 50). When in operation, the outputs from 
the magnetic sensors are processed using the adaptive 
signal processing technique described previously. 
[0077] In known biomagnetic systems, measure- 
ments have to be made in an expensive magnetically 



shielded room. In the arrangement shown in Figure 9, 
the global feedback colls 30a,30b,30c replace this mag- 
netically shielded room. The global feedback coils may 
be configured to feed back the uniform field and also the 
5 first order magnetic field gradients. This provides sub- 
stantial rejection of external field sources. 
[0078] In general, the system may be used for several 
applications where the sensing of magnetic fields is re- 
quired. For example, the system may be particularly 
useful for use in geological surveying, non destructive 
evaluation or in medical applications such as biomag- 
netic scanning. 



1 . A system for measuring one or more magnetic field 
gradient components of a magnetic field compris- 
ing; 

20 

(i) at least two magnetic sensors (25a,25b) for 
sensing a magnetic field, wherein each sensor 
generates a sensor output, said sensor outputs 
having an associated total energy, E, and 

25 wherein at least two of the sensors are ar- 

ranged to sense the magnetic field in substan- 
tially the same direction and 

(ii) means for perfomning adaptive signal 
30 processing of the sensor outputs (31 ,32,30a, 

30b) such that the system is adaptively bal- 
anced, whereby said means generate one or 
more magnetic field gradient components, 



(iii) the system incorporates global feedback 
means (30a, 30b) for providing a substantially 
uniform magnetic field at the two or more mag- 

40 netlc sensors, and 

(iv) the adaptive signal processing means com- 
prise means (34a, 34b, 35) for minimising the 
total energy, E, of the sensor outputs subject to 

45 a constraint, whereby the constraint deter- 

mines which of one or more magnetic field gra- 
dient components is generated. 

2. The system of claim 1 , wherein the means (34a, 
50 34b, 35) for minimising the total energy, E, of the 

sensor outputs may also generate a total magnetic 
field measurement, whereby the constraint deter- 
mines which of one or more magnetic field gradient 
components or a total magnetic field measurement 
55 is generated. 

3. The system of claim 1 , wherein the adaptive adap- 
tive signal processing means include a computer 
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(34a, 34b, 35) having an adaptive signal processing 
algorithm (ASPA) loaded therein. 

4. The system of claim 1 , wherein the global feedback 
means (30a, 30b) are arranged to provide a sub- 
stantially unlfomri magnetic field gradient at the two 
or more magnetic sensors. 

5. The system of claim 1 , wherein the global feedback 
means comprise at least one global feedback coil 

set (30a,30b). 

6. The system of claim 5, wherein the one or more glo- 
bal feedback coil set each comprise two or more 
Helmholtz coils (30a, 30b, 30c). 

7. The system of claim 6, and also including 

means for generating at least one difference 
signal between two sensor outputs, wherein 
said sensor outputs each correspond to a mag- 
netic field in substantially the same direction, 
and 

analogue to digital conversion means (35) for 
converting the one or more difference signal 
and the two or more sensor outputs into equiv- 
alent digital data. 

8. The system of claim 1 , wherein at least one of the 
magnetic sensors is any one of a flux gate, a Hall 
probe or a magneto-resistive sensor 

9. The system of claim 1 , wherein at least one of the 
magnetic sensors is a SQUID magnetometer 

10. The system of claim 9, and also comprising cooling 
means for reducing the temperature of the SQUID 
magnetometer 

1 1 . The system of claim 1 0, wherein each SQUID mag- 
netometer has associated local feedback means 
(26a, 26b) for maintaining a substantially constant 
state of magnetic flux in the respective SQUID mag- 
netometer 

12. The system of claim 1 , comprising; 

at least four magnetic sensors for sensing a 
magnetic field, wherein three of the sensors are 
arranged such that they sense the magnetic 
field in three substantially orthogonal directions 
and 

wherein at least two of said sensors are ar- 
ranged such that they sense the magnetic field 
in substantially the same direction. 

13. The system of claim 12, comprising at least eight 



magnetic sensors. 

14. The system of claim 13, wherein the eight or more 
magnetic sensors are arranged at the vertices of a 

s tetrahedron structure. 

15. The system of claim 14, wherein the global feed- 
back means comprise three global feedback coil 
sets (30a, 30b, 30c) oriented in three substantially 

10 orthogonal directions, for generating a substantially 
unifomi magnetic field at each of the eight or more 
magnetic sensors. 

16. The system of claim 1, comprising at least three 
IS magnetic sensors, wherein at least three of the sen- 
sors are oriented in substantially the same direction 
and whereby the means for performing adaptive 
signal processing may generate a magnetic field 
gradient component of at least second order 

20 

17. The system of claim 1, wherein the system is ar- 
ranged to provide a biomagnetic sensing system 
(55). 

18. The system of claim 17, wherein the global feed- 
back means comprise at least one global feedback 
coil set (30a, 30b or 30c) arranged to sun^ound a 
subject (50) generating a magnetic field to be meas- 
ured, wherein the subject (50) is in close proximity 
with the two or more magnetic sensors (25). 

19. A method for measuring at least one magnetic field 
gradient component using a gradiometer compris- 
ing the steps of; 

(1) sensing a magnetic field component at two 
or more positions using two or more magnetic 
sensors (25a, 25b), wherein at least two of the 
magnetic field components are sensed in sub- 
stantially the same direction, 

(ii) generating two or more output signals, hav- 
ing an associated total energy, E, correspond- 
ing to said magnetic field components, 

(ill) providing global feedback in the form of a 
substantially unifomn magnetic field at the two 
or more magnetic sensors (25a, 25b), 

(Iv) performing adaptive signal processing of 
the output signals and minimising the total en- 
ergy, E, of the output signals subject to a con- 
straint such that the gradiometer is adaptively 
balanced, 

(v) constraining the minimisation of the total en- 
ergy, E, such that a magnetic field gradient 
component is generated, and 
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(vi) generating at least one magnetic field gra- 
dient component measurement. 

20. The method of claim 1 9, comprising the steps of; 

5 

(i) minimising the total energy, E, of the output 
signals subject to a constraint such that the 
gradiometer is adaptively balanced, and 

(ii) constraining the minimisation of the total en- io 
ergy, E, such that a total magnetic field meas- 
urement Is generated. 

21 . The method of claim 1 9, wherein the magnetic field 
components at two or more positions are sensed is 
using two or more magnetic sensors (25a, 25b), and 
further comprising the step of providing a substan- 
tially uniform magnetic field gradient at the two or 
more magnetic sensors. 

20 

22. The method of claim 21 for measuring magnetic 
field components in a human subject (50) in close 
proximity with the two or more magnetic sensors 
(25). 

25 

23. The method of claim 19, comprising the steps of; 

(i) sensing the magnetic field component at 
three or more positions, wherein at least three 

of said magnetic field components are sensed 30 
in substantially the same direction, 

(ii) minimising the total energy, E, of the output 
signals subject to a constraint such that the 
gradiometer is adaptively balanced, and 35 

(ill) constraining the minimisation of the total 
energy, E, such that a magnetic field gradient 
component of at least second order is generat- 
ed. 40 



PatentansprQche 

1 . System zum Messen einer oder mehrerer Magnet- 
feldgradientenkomponenten eines Magnetfeldes 
mit: 

(i) wenigstens zwei Magnetsensoren (25a, 
25b) zum Erfassen eines Magnetfeldes, wobei so 
jeder Sensor einen Sensorausgang erzeugt, 

die Sensorausgange eine dazugehdrige Ge- 
samtenergie E aufweisen und wenigstens zwei 
der Sensoren angeordnet sind, um das Ma- 
gnetfeld in im wesentlichen derselben Richtung ss 
zu erfassen, und 

(ii) einer Vorrichtung zum Durchfuhren der ad- 
aptiven Signalverarbeitung der Sensorausgan- 



ge (31, 32, 30a, 30b), so daB das System ad- 
aptiv ausbalanciert ist, wobei die Vonrichtung 
eine oder mehrere Magnetfeldgradientenkom- 
ponenten erzeugt, 

dadurch gekennzeichnet, daB 

(iii) das System eine globale Riickkopplungs- 
vorrichtung (30a, 30b) zum Erzeugen eines im 
wesentlichen gteichfonnigen Magnetfeldes bei 
den zwei oder mehreren Magnetsensoren um- 
faf3t und 

(iv) die adaptive Signalverariseitungsvorrich- 
tung Vorrichtungen (34a, 34b, 35) zum Mini- 
mieren der Gesamtenergie E der Sensoraus- 
gange in Abhangigkeit von einer Einschran- 
kung umfaBt, wobei durch die Einschrankung 
bestimmt wird, welche der einen oder mehre- 
ren Magnetfeldgradientenkomponenten er- 
zeugt wird. 

2. System nach Anspruch 1, bei dem die Vorrichtun- 
gen (34a, 34b, 35) zum Minimieren der Gesamten- 
ergie E der Sensorausgange auBerdem eine Ge- 
samtmagnetfeldmessung erzeugen, wobei durch 
die Einschrankung bestimmt wird, ob eine oder 
mehrere Magnetfetdgradientenkomponenten er- 
zeugt werden Oder ob eine Gesamtmagnetfeldmes- 
sung durchgefuhrt wird. 

3. System nach Anspruch 1 , bei dem die adaptive Si- 
gnalverariseitungsvorrichtung einen Computer 
(34a, 34b, 35) nebst einem adaptiven Slgnalverar- 
beitungsalgorithmus (ASPA) darin umfaBt. 

4. System nach Anspruch 1. bei dem die globalen 
Ruckkopplungsvon'ichtungen (30a, 30b) angeord- 
net sind, um einen im wesentiichen gleichfonmigen 
Magnetfetdgradienten an den zwei oder mehr Ma- 
gnetsensoren zu erzeugen. 

5. System nach Anspruch 1, bei dem die globale 
Ruckkopplungsvorrichtung wenigstens einen glo- 
balen Ruckkopplungsspulensatz (30a, 30b) um- 
faBt. 

6. System nach Anspruch 5 , bei dem der eine oder die 
mehreren globalen Ruckkopplungsspulensatze je- 
weils zwei oder mehr Helmholtz-Spulen (30a, 30b. 
30c) umfassen. 

7. System nach Anspmch 6. das auBerdem umfaBt: 

eine Vorrichtung zum Erzeugen wenigstens ei- 
nes Differenzsignals zwischen zwei Sensor- 
ausgangen, wobei die Sensorausgange jeweils 
einem Magnetfeld in im wesentlichen der gtei- 
chen Rbhtung entsprechen, und 
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eine Analog/Digital-Wandlervorrichtung (35) 
zumWandeIn deseinen oderdermehreren Dlf- 
ferenzsignale und der zwei oder mehreren 
Sensorausgange in aquivalente Digitaldaten. 

8. System nach Anspruch 1 , bel dem wenigstens einer 
der Magnetsensoren ein MagnetfluB-, ein Hall-Sen- 
sor Oder ein Magnetwiderstandssensor 1st. 

9. System nach Anspruch 1 , bei dem wenigstens einer 
der Magnetsensoren ein SQUID- Magnetometer ist. 

10. System nach Anspruch 9, das auBerdem eine Kiihl- 
vorrichtung zum Absenken der Temperatur des 
SQUID-Magnetometers umfaBt. 

11. System nach Anspruch 10, bei dem jedes SQL) ID- 
Magnetometer eine dazugehorige loltale Riick- 
kopplungsvorrichtung (26a, 26b) fiir das Aufrecht- 
erhalten eines im wesentiichen konstanten Zustan- 
des des Magnetflusses In dem Jeweiilgen SQUID- 
Magnetometer umfaBt. 

12. System nach Anspruch 1 mit wenigstens vier Ma- 
gnetsensoren fiir das Erfassen eines Magnetfel- 
des, wobei drei der Sensoren derart angeordnet 
sind, daB sie das Magnetfeld in drei im wesentii- 
chen orthogonalen Richtungen erfassen und wobei 
wenigstens zwei der Sensoren so angeordnet sind, 
daB sie das Magnetfeld in im wesentiichen der glei- 
Chen RIchtung erfassen. 

13. System nach Anspruch 1 2 mit wenigstens acht Ma- 
gnetsensoren. 

14. System nach Anspruch 13, bei dem die acht oder 
mehr Magnetsensoren an den Spitzen einer Tetra- 
ederstruktur angeordnet sind. 

15. System nach Anspruch 14, bei dem die globale 
Ruckkopplungsvorrichtung drei globale Ruckkopp- 
lungsspulensatze (30a, 30b, 30c) umfaBt, die in drei 
im wesentiichen orthogonale Richtungen orientiert 
sind, urn ein im wesentiichen glelchfonniges Ma- 
gnetfeld an jedem der acht oder mehr Magnetsen- 
soren zu erzeugen. 

16. System nach Anspruch 1 mit wenigstens drei Ma- 
gnetsensoren, bei dem wenigstens drei der Senso- 
ren in im wesentiichen der gleichen RIchtung orien- 
tiert sind und bei dem die Vorrichtung fur die Durch- 
fuhrung der adaptiven SIgnalverarbeitung eine Ma- 
gnetfeldgradientenkomponente von wenigstens 
zweiter Ordnung erzeugt. 

17. System nach Anspruch 1, bel dem das System so 
angeordnet ist, daB es ein Biomagnetfeldsensorsy- 
stem (55) bildet. 



18. System nach Anspruch 17, bei dem die globale 
Ruckkopplungsvorrichtung wenigstens einen glo- 
balen Rucl<kopplungsspulensatz (30a, 30b oder 
30c) umfaBt, der so angeordnet ist, daB er ein Ob- 
5 jekt (50) umgibt. das ein zu messendes Magnetfeld 
erzeugt, wobei das Objekt (50) sich in unmittelbarer 
Ndhezu den zwei oder mehreren Magnetsensoren 
(25) befindet. 

10 19. Verfahren zum Messen wenigstens einer Magnet- 
feldgradientenkomponente unter Venwendung ei- 
nes Gradlometers mit den Schrltten: 
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(i) Erfassen einer Magnetfeldkomponente an 
zwei Oder mehreren Positionen unter Verwen- 
dung von zwei oder mehreren Magnetsensoren 
(25a, 25b) , wobei wenigstens zwei der Magnet- 
feldkomponenten in Im wesentiichen der glei- 
chen Rtehtung erfaBt werden, 

(ii) Erzeugen von zwei oder mehreren Aus- 
gangsslgnalen mit einer dazugehorigen Ge- 
samtenergie E, die den Magnetfeldkomponen- 
ten entspricht, 

(iii) Erzeugen einer globalen Ruckkopplung in 
Form eines im wesentiichen gleichfomiigen 
Magnetfeldes an den zwei oder mehreren Ma- 
gnetsensoren (25a, 25b), 

(iv) Durchfuhren der adaptiven SIgnalverarbei- 
tung der Ausgangssignale und Minimieren der 
Gesamtenergie E der Ausgangssignale in Ab- 
hangigkelt von einer Einschrankung, so daB 
das G radiometer adaptiv ausbalanclert ist, 

(v) Einschranken der Minlmterung der Gesam- 
tenergie E, so daB eine Magnetfeidgradienten- 
komponente erzeugt wird, und 

(vi) Erzeugen wenigstens einer Magnetfeldgra- 
dientenkomponentenmessung. 

20. Verfahren nach Anspruch 19 mit den Schritten: 

(i) Minimieren der Gesamtenergie E der Aus- 
gangssignale in Abhangigkeit von einer Ein- 
schrankung, so daB das Gradiometer adaptiv 
ausbalanciert ist, und 

(ii) Einschranken der MInimierung der Gesam- 
tenergie E, so daB eine Gesamtmagnetfeld- 
messung erzeugt wird. 



21. Verfahren nach Anspruch 19, bei dem die Magnet- 
ic feldkomponenten an zwei oder mehreren Positio- 
nen unterVerwendung von zwei oder mehreren Ma- 
gnetsensoren (25a, 25b) erfaBt werden und das au- 
Berdem als Schritt das Erzeugen eines im wesent- 
iichen gleichfomiigen Magnetfeldgradientenan den 
zwei Oder mehreren Magnetsensoren umfaBt. 

22. Verfahren nach Anspruch 21 zum Messen von Ma- 
gnetfeldkomponenten in einem menschlichen Ob- 
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jekt (50) in unmittelbarer Nahe zu den zwel Oder 
mehreren Magnetsensoren (25). 

23. Verfahren nach Anspruch 19 mit den Schritten: 

(i) Erfassen der Magnetfeldkomponentean drei 
Oder mehreren Positionen, wobei wenigstens 
drel derMagnetfeidkomponenten in im wesent- 
lichen dergletehen RIchtung erfa3t werden, 
(ti) Minimleren der Gesamtenergie E der Aus- 
gangssignale in Abhangigkeit von einer Ein- 
schrankung, so daB das Gradionneter adaptiv 
ausbalanciert wird, und 
(ill) Einschranken der Minimiemng der Gesam- 
tenergie E, so da3 eine Magnetfeldgradienten- 
komponenten von wenigstens zweiter Ordnung 
erzeugt wird. 



Revendications 

1 . Syst^me pour mesurerune ou plusieurs composan- 
tes de gradient de champ magn^tique d'un champ 
magndtlque comportant : 

(i) au moins deux capteurs magndtiques (25a, 
25b) pour detecter un champ magn^tique, cha- 
que capteur generant une sortie de capteur, 
lesdites sorties de capteur ayant une energie 
totale associ6e, E, et au moins deux des cap- 
teurs 6tant census pour detecter le champ ma- 
gn^tique pratiquement dans la m&me direction, 
et 

(ii) des moyens pour effectuer un traitement de 



g6n6rer une mesure de champ magn^tique total de 
sorte que la contrainte d^tennine qui pamfii une ou 
plusieurs composantes de gradient de champ ma- 
gnetique ou une mesure de champ magnetique to- 
ff tal est gSnerde. 

3. Syst&me selon la revendteation 1 , dans lequel les 
moyens de traitement de signal adaptatif compor- 
tent un ordinateur (34a, 34b, 35) ayant un algorith- 

10 me de traitement de signal adaptatif (ASPA) charge 
a I'interieur. 

4. Syst^me selon la revendicatlon 1 , dans lequel les 
moyens de reaction globale (30a, 30b) sont congus 

15 pour foumir un gradient de champ magnetique pra- 
tiquement uniforme au niveau des deux ou plus de 
deux capteurs magn^tiques. 

5. Systdme selon la revendicatlon 1 , dans lequel les 
20 moyens de rdaction globate comporte au moins un 

ensemble de bobines de reaction globale (30a, 
30b). 

6. Syst^me selon la revendicatlon 5, dans lequei le ou 
25 les ensembles de bobines de reaction globale com- 

portent chacun deux ou plus de deux bobines de 
Hefmholtz (30a, 30b, 30c). 

7. Le systdme selon la revendicatlon 6, et comportant 
30 6galement 

des moyens pour g6n6rer au moins un signal 
diffdrentiet entre deux sorties de capteur, lesdi- 
tes sorties de capteur correspondant chacune 
^ un champ magnetique pratiquement dans la 
meme direction, et 

des moyens de conversion analogique-nume- 
rlque (35) pour convertir le ou les plusieurs si- 
gnaux diff^rentiels et les deux ou plus de deux 
sorties de capteur en donn^es num^riques 
dquivalentes. 



signal adaptatif des sorties de capteur (31 , 32, 35 
30a, 30b) de sorte que le systeme est equilibre 
de mani^re adaptative, de sorte que lesdits 
moyens gdn^rent une ou plusieurs composan- 
tes de gradient de champ magnetique, 

40 

caracterise en ce que 



(iii) le systeme comporte des moyens de reac- 
tion globale (30a, 30b) pour foumir un champ 
magnetique pratiquement unifomie au niveau 45 
des deux ou plus de deux capteurs magneti- 

ques, et 

(iv) les moyens de traitement de signal adapta- 
tif comportent des moyens (34a, 34b, 35) pour 
minimiser renergie totale, E, des sorties de so 
capteur soumises ^ une contrainte, de sorte 
que la contrainte determine quelle composante 
parmi une ou plusieurs composantes de gra- 
dient de champ magnetique est generee. 

55 

2. Systeme selon la revendicatlon 1 , dans lequel les 
moyens (34a, 34b, 35) pour minimiser i'energie to- 
tale, E, des sorties de capteur peuvent egalement 



8. Systeme selon la revendicatlon 1 , dans lequel au 
moins un des capteurs magnetiques est un dispo- 
sitif quelconque parmi un detecteur de flux, une 
sonde de Hall ou un capteur magneto-resistif. 

9. Systeme selon la revendicatlon 1 , dans lequel au 
moins un des capteurs magnetiques est un magne- 
tometre SQUID. 

10. Systeme selon la revendicatlon 9, et comportant 
egalement des moyens de refroidissement pour re- 
duire la temperature du magnetometre SQUID. 

11. Systeme selon la revendicatlon 10. dans lequel 
chaque magnetometre SQUID a des moyens de 
reaction locale associes (26a, 26b) pour maintenir 
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un etat pratiquement constant de flux magn^tique 
dans le magnetonndtre SQUID respectif. 

12. Syst^me selon la revendication 1 , comportant : 

au nnoins quatre capteurs magn6tiques pour s 
d^tecter un champ magn^tique, trols des capteurs 
6tant conpus de sorte qu'lfs d^tectent le champ ma- 
gnetique dans trols directions pratiquement ortho- 
gonales et au moins deux desdits capteurs sont 
congus de sorte qu'Hs detectent le champ magn6ti- io 
que pratiquement dans la meme direction. 

13. Systeme selon la revendication 12, comportant au 
moins huit capteurs magn^tiques. 



14. Systeme selon la revendication 13, dans lequel les 
huit ou plus de huit capteurs magnetiques sont 
agences au niveau des sommets d'une structure de 
t6tra&dre. 
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15. Systdme selon la revendication 14, dans lequel les 
moyens de reaction globale comportent trols en- 
sembles de bobines de reaction globale (30a, 30b, 
30c) orientes selon trols directions pratiquement or- 
thogonales, pour generer un champ magnetique 25 
pratiquement unifomne au niveau de chacun des 
huit ou plus de huit capteurs magnetiques. 

16. Systeme selon la revendication 1, comportant au 
moins trois capteurs magnetiques, dans lequel au 20 
moins trols des capteurs sont orlent6s pratique- 
ment dans la meme direction et de sorte que les 
moyens pour effectuer un traitement de signal 
adaptatif peuvent g^n^rer une composante de gra- 
dient de champ magnetique au moins du second ss 
ordre. 

17. Systeme selon la revendication 1, dans lequel le 
systeme est con^u pour foumir un systdme de de- 
tection blomagn^tique (55). 40 

18. Systeme selon la revendication 17, dans lequel les 
moyens de reaction globale comportent au moins 
un ensemble de bobines de reaction globale (30a, 
30b ou 30c) agence pour entourer un sujet (50) ge- 4s 
nerant un champ magnetique k mesurer, dans le* 
quel le sujet (50) est k proximity 6troite des deux 

ou plus de deux capteurs magnetiques (25). 

19. Procede pour mesurer au moins une composante so 
de gradient de champ magnetique utilisant un gra- 
diomanometre, comportant les etapes consistant 

k: 

(i) detecter une composante de champ magne- ss 
tique au niveau de deux ou plus de deux posi- 
tions en utilisant deux ou plus de deux capteurs 
magnetiques (25a, 25b), au moins deux des 



composantes de champ magnetique etant de- 
tectees pratiquement dans la meme direction, 
(il) generer deux ou plus de deux signaux de 
sortie, ayant une energie totale associee, E, 
correspondant auxdites composantes de 
champ magnetique, 

(Hi) foumir une reaction globale ayant la fomne 
d'un champ magnetique pratiquement unifor- 
me au niveau des deux ou plus de deux cap- 
teurs magnetiques (25a, 25b), 

(iv) effectuer un traitement de signal adaptatif 
des signaux de sortie et minlmiser Tenergie to- 
tale, E, des signaux de sortie soumis a une con- 
trainte de sorte que le gradiomanometre est 
equillbre de maniere adaptatlve, 

(v) imposer la minimisation de i'energie totale, 
E, de sorte qu'une composante de gradient de 
champ magnetique est generee, et 

(vi) generer au moins une mesure de compo- 
sante de gradient de champ magnetique. 

20. Precede selon la revendication 1 9, comportant les 
etapes consistant k : 

(i) minlmiser I'energie totale, E, des signaux de 
sortie soumis a une contrainte de sorte que le 
gradiomanometre est equilibre de maniere 
adaptative, et 

(II) imposer la minimisation de I'energie totale, 
E, de sorte qu'une mesure de champ magneti- 
que total est generee. 

21. Precede selon la revendication 19, dans leque! les 
composantes de champ magnetique au niveau de 
deux ou plus de deux positions sont detectees en 
utilisant deux ou plus de deux capteurs magneti- 
ques (25a, 25b), et comportant de plus I'etape con- 
sistant k fournir un gradient de champ magnetique 
pratiquement unifonne au niveau des deux ou plus 
de deux capteurs. 

22. Precede selon la revendication 21 , pour mesurer 
des composantes de champ magnetique dans un 
sujet humain (50) a proximite etroite des deux ou 
plus de deux capteurs magnetiques (25). 

23. Precede selon la revendication 19, comportant les 
etapes consistant k : 

(I) detecter la composante de champ magneti- 
que au niveau de trois ou plus de trois posi- 
tions, au moins trois desdites composantes de 
champ magnetique etant detectees pratique- 
ment dans la m§me direction, 

(II) minlmiser Tenergie totale, E, des signaux 
des soumis k une contrainte de sorte que le 
gradiomanometre est equilibre de maniere 
adaptative, et 
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(iii) imposer la minimisation de I'energie totale, 
E, de sorte qu'une composante de gradient de 
champ magnetique au moins du second ordre 
est generee. 
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